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We study the configurational structure of the point-island model for epitaxial growth in one 
dimension. In particular, we calculate the island gap and capture zone distributions. Our model is 
based on an approximate description of nucleation inside the gaps. Nucleation is described by the 
joint probability density p^^ {x,y), which represents the probability density to have nucleation at 
position X within a gap of size y. Our proposed functional form for p^^ {x^y) describes excellently 
the statistical behavior of the system. We compare our analytical model with extensive numerical 
simulations. Our model retains the most relevant physical properties of the system. 
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I. INTRODUCTION 

The study of the theory of the submonolayer film 
growth is of great general interest, in part because its 
application in the construction of microelectronic de- 
vices. In epitaxial growth, monomers are deposited onto 
a substrate at a constant rate. The monomers diffuse 
until they are captured by an island or another diffu- 
sive monomer. The islands are just clusters of immobile 
monomers. The size of an island depends on the number 
of monomers which have attached to it. Thus, the is- 
land size increases in time due to the capture of diffusing 
monomers. These kinds of systems exhibit many inter- 
esting non-equilibrium phenomena, as discussed in, e.g., 
Refs. [Il[2] 

There are basically two kinds of models which have 
been developed to study these systems: In the extended 
island model, the islands occupy more than one site on 
the lattice, and their shape is not trivial |2H7j. In the 
simpler point-island model, the islands just occupied one 
site in the lattice, and their size is simply the number 
of monomers which belong to the island P4T6]. In this 
paper, we focus on the point-island model, which is very 
accurate for low coverages. A major advantage of this 
simplification is that it produces better statistics than 
the extended island model. 

In our model only the monomers are mobile, the is- 
lands formed by two or more monomers are completely 
static (stable). In the literature this condition is usu- 
ally denoted "critical nucleus size" i = 1 (where i is the 
size of the largest unstable island). In our irreversible- 
growth model, monomers must land or hop onto an al- 
ready occupied site in order to be incorporated |8l [IS] ; in 



an alternative model features incorporation of monomers 
arriving at sites adjacent to the island [16 . 

The quantities that are commonly used to describe the 
evolution of the point-island model are the density of 
monomers, A^i; the density of islands with size j, Nj] 
the rate of deposition of monomers, F; and the diffusion 
rate of monomers D. The evolution of this system is fre- 
quently described in terms of the coverage = Ft^ where 
t is the time. We restrict our studies to the aggregation 
regime, where there is a quasi steady-state. 

We are particularly interested in the spacing (gap) dis- 
tribution functions between islands p^^\S)^ the capture 
zone distribution P{S)^ and the island-island pair corre- 
lation function G(r). As usual, p''^\S)dS is the probabil- 
ity that for an island at the origin we find another island 
at a distance between S and S + dS^ with the condition 
that there are k additional islands inside the gap between 
them. The standard definition for the scaled spacing is 
s = S/ (5), with (S) the average of S. The scaled spacing 
distributions are given by 



p('=)(s) = (5)pW(5(5)) 



(1) 



On the other hand, we use the definition for P{S) given 
in [8 . Then, in our one-dimensional (ID) system, the 
capture zone of an island is simply the distance between 
the midpoints of the gaps to the left and to the right of 
the island. 

The functional forms of p^^\s) and P{s) for arbitrary 
dimension and critical nucleus size have been the subject 
of recent discussion and some controversy [HI [TTHTQ] . A 
particular issue is whether the generalized Wigner sur- 
mise (GWS) 
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adequately describes the distribution and the reliability 
the simple relationship between j3 and i deduced with 
mean field [TT^ and later refined with more sophisticated 
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arguments [19] . There is still no general consensus about 
them even in the case of i = 1 in one dimension (ID). 
In this paper we calculate analytically and numerically 
the spacing distribution functions p^^\s) for /c > 0, a 
meaningful concept for capture zones in ID but less so 
in higher dimensions. All previous studies focus just in 
p^^\s) [ll[T0l|T7]. We specially focus in their functional 
forms in the limit of large and small values of s, where 
deviations from the form of Eq. Q were observed in 
painstaking computations in ID and 2D, see Ref. [8]. We 
propose an analytical model to find approximate expres- 
sions for these functions. Our model is based on a de- 
tailed description of the nucleation mechanism, following 
principally Blackman and Mulheran [T], hereafter BM. 
We also calculate the pair correlation function G(r). In 
all cases our model is compared with several numerical 
simulations. This paper is organized as follows: In Sec. 2 
we discuss briefiy some implications of the ID model. In 
Sec. 3 we provide an accurate description of the nucle- 
ation process. In Sec. 4 and Sec. 5, we calculate approxi- 
mately the spacing distribution functions, while in Sec. 6 
we calculate the island-island pair correlation function. 
In Sec. 7, we assess the viability of applying Eq. ([2| to 
experimental data. In Sec. 8, we give conclusions. 



From its definition, it is clear that p'^^^si +52) is related 
to the capture zone distribution P{{si-\- S2) /2) according 
to 



P{s) = dxp^^\x)8(s - I) = 2p^'H25), (4) 

where x = 5i + 52. 

Thus, the capture zone distribution and the next- 
nearest spacing distributions are equivalent in one di- 
mension. 



Si 



FIG. 1. Relation between P{s) and p'^^\s). Black blocks are 
islands, while white regions are gaps. 



II. ID POINT ISLAND MODEL 

In the ID case, we have a ring divided into independent 
sections called gaps. Each gap starts and ends with an 
island. Within the gaps, there may be several monomers 
performing random walks. A monomer inside of a par- 
ticular gap must eventually either merge with one of the 
islands at the ends of the gap or combine with another 
monomer to nucleate a new island. In no case can the 
monomer reach a different gap. 

The spatial distribution of the islands, i.e., the distri- 
bution of the sizes of the gaps, is given by the spacing 
distribution functions p^^\s). The simplest case corre- 
sponds to the nearest-neighbor distribution p*^^^ (s), which 
represents the probability density to find a gap with an 
scaled size s between two islands with no additional is- 
land in between. 

The pair correlation function G(r) is related to p^^\s) 

by 

00 

G(r) = ^pW(r). (3) 

k=Q 

As noted, the length of the capture zone of an island 
is simply the distance between the midpoints of the gaps 
to the left and to the right of the island. In our ID 
model, the next-nearest neighbor distribution p^^\s) is 
related to P{s) as follows. Consider Fig [l] where the 
black squares represent islands. Let Si + 5*2 be the dis- 
tance between the islands A and B. Then p^^\Si + ^'2) 
is the probability density to find a gap with size + ^2 
given that there is an additional island inside the gap. 



III. SPATIAL DESCRIPTION OF THE 
NUCLEATION 

In order to find an expression for p^^\s)^ we must de- 
scribe the creation mechanism of new gaps. To reach an 
appropriate description of the nucleation, we make the 
following definitions: Let y) be the joint probabil- 

ity density that a given new nucleation occurs at position 
X inside a gap of length y. Of course x < y; otherwise 
Pn^{x^ y) = 0. Let p^{x) be the probability density that 
a given nucleation occurs at position x inside a gap of any 
size. Then 

p^ix)= / dyp^''ix,y). (5) 

J X 

Similarly, the probability density p^ (y) that a given nu- 
cleation occurs anywhere inside a gap of length y is given 
by 

pl{y)= r dxp^^'ix.y). (6) 
Jo 

Finally, the probability density that, a nucleation occurs 
at position x in a gap of length y under the condition 
X = x/y is 

pt{\) = dy £ dxp^ix, y)5(^X- ^) 

POO 

= / dyyp^''iXy,y). (7) 

The formation of a new island depends directly on the 
probability that two monomers reach the same site of 
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FIG. 2. (Color online) Behavior of (a) the normalized reduced density n(A = x/y) (see Eq. (10)) and (b) the probability 
p^{X = x/y)) that a nucleation occurs at position x in a gap of length y, for different values of R. The agreement between the 
analytical model and the numerical data is excellent overall. There are some slight differences near the maximum of p^{X). 



the lattice at the same time. This quantity is of course 
related to the density of monomers inside the gap. In 
BM, the average density of monomers, ni(x,?/), inside of 
a gap of length y was approximated by using its expres- 
sion in the stationary state. In this regime there is quasi- 
equilibrium between the deposition of new monomers and 
their merging with existing islands. We neglect the de- 
position onto the top of islands. Then ni{x^y) satisfies 
a one-dimensional diffusion equation with the boundary 
conditions ni(0,?/) = ni{y^y) = 0. This leads to 



ni{x,y) = {2R) ^x{y - x), 



(8) 



with R = D/F [1 . By using this expression and the 
assumption that the probability of a new nucleation at x 
is proportional to ni(x,?/)^, we have 



(9) 



where A = x/y [1 . While we have neglected the inter- 
action among monomers within a gap, a more refined 
estimate of p^(A) can be made by taking them into ac- 
count and using the full multiparticle density [20l|2T]. We 
tested numerically the validity of this approximation. As 
seen in Fig. |2j there is excellent agreement between the 
numerical data and the analytical equations. Note that 
in Fig. [2] we plotted the normalized reduced density n(A) 
given by 



n(A) 



poo ny 
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(10) 

instead of ni{x^y). Here /is is the 3^^ moment of p^^\y) 
and the factor 12 R/ fi^, guarantees the correct normaliza- 
tion of n(A) in the interval < A < 1. From its definition, 
it is clear that n{\) gives the average density of monomers 
inside of a gap at the relative position \ = x/y. A major 



advantage of n(A) is that it does not depends explicitly 
on the length of the gap, in contrast to ni(x, y). The nu- 
merical simulation of the point-island model was carried 
out following the standard procedures [T| [9l [10] . In our 
simulations we took a lattice with length L = 2 x 10^ 
sites (with periodic boundary conditions) and gathered 
statistics from over 5000 realizations. 

Since p^(A) vanishes for A = and 1, a new island 
is unlikely to form near the boundaries of a gap. This is 
the origin of the effective repulsive force between adjacent 
islands. As we shall shortly see, this implies that for small 
values of 5, p^^\s) ex . We expect that a = 2 because 
p^(A)(xA2for A^O. 

BM also proposed that p^{y) / Fm^v) oc y^ , where 
Fuiy) is the number of gaps with size y given that there 
are M gaps. Fm{S) is naturally related to p^^\s) by 
Fm{S) = {W^ /L)p^^\s). We calculate numerically this 
quotient for different values of R. As shown in Fig. [3) our 
numerical results suggest that p^{y)/FM{y) oc y^ with 
7 ^ 3 for 5 > 1.7 and 7 ^ 4 for 5 < 1.7 rather than 7 = 5. 
A similar result was found in Ref. [16] for submonolayer 
deposition in 2D. Their numerical results suggests that 
7(5) ^ (4 + 5)/(2 + s). Hence, for small gaps we have 
7^2 while for large gaps 7^1. 

Implicitly in BM, the probability density p^{y) was 
written as 



dxnl{x,y)]p^^\y) ex y^p^^\s). (11) 



Underlying Eq. (11) are the following two approxima- 
tions: First, the integral is based on the law of mass 
action; Fig. [2] justifies this. Second, BM supposed that 
p^ {y) can be written as the product of the probability to 
have a nucleation inside a gap of size y and the number of 
gaps FM{y) with this size. Our numerical results, shown 
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in Fig. [3j do not support that simplification. 



IV. NEAREST-NEIGHBOR SPACING 
DISTRIBUTION p^^\s) 
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FIG. 3. (Color online) Behavior of the quotient Pn (s) / P^'^\s) . 
It seems that there are two different regimes, one for small 
gaps and the other for large gaps. 

Following our previous observations, Pn^{x^y) can be 
written as the independent product of the probability 

{y) to nucleate inside a gap with size y and the condi- 
tional probability Pn{x/y)/y that this nucleation occurs 
at X. After some algebra we find 



30 

^7 



(12) 



where ji^ is the 7^^ moment of p^^\y). It is easy to show 
that Eq. ([l2| satisfies the form of p^(A) given in BM. The 
condition oc y'^ is satisfied as well. From 

Eqs. ([5|, ([6| and ([l2|, it is straightforward to show 



pI{s) 



(13) 



and 
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1: 



y 



dy 

5-7 



p^'\y){y-sf. (14) 



Equations (13) and (14) describe the mechanism of for- 
mation of new islands. Equation (13) shows that the 



probability of nucleation inside a gap of size y is the prod- 
uct of the number of such gaps and y^ . As we will see 



7 7^ 5 in contradiction to Eq. (11). Equation (14) shows 
how Pn (s) depends on the number of monomers inside 
the gap; for small s, p^(s) oc 5^, consistent with the law 
of mass action. The vanishing of {y — s)'^ diS y ^ s and 
as s ^ reflects that nucleation does not occur at the 
gap borders because ipso facto the monomer concentra- 
tion vanishes there. 



Following BM, it is possible to find an analytical ex- 
pression which relates p^^\y) with p^i^) and p^(s): The 
effect of one single nucleation is described by 



Fm+i{S) - Fm{S) = -pI{S) + Rn{S), (15) 



with 



Rn{S) = dy dxp^'^ix, y) {5{x - S) ^ 5{x - y ^ S)) 

Jx Jo 

= / dy{p^^{S,y)+p^^{S-y,y)). (16) 

Due to the symmetry of p^{X) around its maximum, 
pr(S,y)=pr{S-y,y). Thus, 



FM+iiS) - Fm(5) = -p'^iS) + 2p^{S). 



(17) 



Finally, taking into account the relation between Fm{S) 
and p*^^^(s), we find 

s*^ + 2p(°)(«) = -pr(«) + 2p;f(.). (18) 



From Eqs. (13), (14) and (18) we can write 
dp^^'^s) 



ds 



2^^)p^^\s) = 

^7 / /^7 



/ 

J s 



{y-sf 



dy^^^/'Hy)- 



(19) 

This integro-differential equation can be written in the 
following integral form 



p(°)(s) = 



2 e 



f 

Jo 



dyyPn{y)e' 



(20) 



where p^iv) is given by Eq. (14). BM attributed the 
exponential tail of p'^^^s) to the fact that p^is) cx y^ , 
regardless of the form of p^(A). We recall that, in their 
case, 7 = 5. Consequently, they claim that, for large 
values of s, p'^^^s) cx exp(— 5^/5 /i^). Furthermore, the 
exponential tail of p*^^^ (y) depends on the value of 7 in 
the T atio p^{s)/p^^\s) cx regardless the form of p^(A). 
This means that the fragmentation process for large val- 
ues of s depends on the probability of choosing the gap 
to fragment, while for small values of 5, it depends ex- 
clusively on the probability of choosing the place of frag- 
mentation. 

From Eqs. ([l4|) and ( [2Q| ) it can be shown that in the 
limit s <^ 1, p^{s) = (T5//i^)s2 + 0{s^) and p^ (s) = 
(30//i^)s^ + O(s^). In the opposite limit 5 ^ 00 we have 
p^^\s) cx exp(— s^/7/i^) (see Appendix). 

As noted in BM, an integro-differential equation like 
Eq. ([l9| is hard to solve analytically or numerically. 
However, it can be solved easily in an statistical way: 
A line of length L is taken and an array of points on 
the line is generated. Each new point is introduced via 



5 




10 15 20 25 




10-^ 
10-^ 
10-^ 
10-^ 



■ — Model y' 

■ - - ■ Model y"* 

Model y^ 

TRM 

. R=10', 6=0.5 
• R=0.5x10', 9=1 
. R=10', 9=0.1 



(a) 



10 15 20 



10^ 

10-^ 
10-^ 
10-^ 
10-^ 
10-^ 



25 



FIG. 4. (Color online) Comparison of the statistical solution of Eq. (20) and the numerical simulation of the point-island model 
with i = 1. In (a) we show the first four spacing distribution functions and (b) shows the behavior of p^'^^s) for large values of 
s. 
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FIG. 5. (Color online) The behavior of p^^\s) and (s) are shown in (a) and (b), respectively. There is excellent agreement 
between Eqs. (24) and (25), respectively, and the numerical data from the simulation of the point-island model. 



weighted random number generation. First, a pair of ex- 
isting points is selected between which to introduce the 
new point. The selection is weighted by the 7^^ power 
of the separation of the points. We used 7 = 3,4 and 
5. Then, the actual position in the gap for insertion x 
is selected with weighting given by the function p^{x/y). 
We use L = 10000 and gather data for 100 particles and 
50000 realizations. Henceforth, we call this procedure 
the statistical solution of Eq. (19). 



We also seek a simple approximate expression for 
p'^^^s) which can be used instead of the statistical so- 
lution of the integro-differential equation. In order to 
find it, we focus on the known moments of p^^\s). The 
first two moments are chosen to satisfy the standard nor- 



malization conditions. In Ref. [T], it was shown that in 
the aggregation regime the density of monomers, A^i, is 
related to the islands density, TV, according to 



UN^NiR^ 1.6, 



(21) 



where /is is the third moment of p^^\s). There is an- 
other condition which can be extracted from p^^^(s). In 
the aggregation regime, the spacing distribution func- 
tions have only a weak dependence on time. Because 
of this we can interpret the point-island model as a one- 
dimensional system in which the particles (islands) inter- 
act with their nearest neighbors under a potential v{s). 
In the limit of small coverage and small values of 5, we 
can write p'^^^s) ~ exp[—v{s)/kBT]. For simplicity, we 
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henceforth set the temperature to he ksT = 1. It fol- 
lows that v{s) ^ —2 ln{s). It is easy to check the validity 
of this approximate equation from our results for p^^\s). 
The average interaction energy per particle u is constant; 
using the maximum entropy method (MEM) [22-.. - -24^ , we 
can approximate u by 



1 

.(0) 



dS v{S) p[J^^ 



(S) 



(22) 



To determine PuEui^)^ tdike note of four independent 
conditions: two normalization conditions, the third mo- 
ment of p^^\s) and the average of v{s). Then PuEui^) 
is given by 



^Ci+C2S+C3S^+C4ln(s2) 



(23) 



where Ci are constants. We conclude that an approxi- 
mate analytical form for p^^^ (s) is 



with A, B and C constants to be determined. 



(24) 



finding straightforwardly 



By using Eqs. (18) and (24), we can calculate p^{s)^ 



ff(°)(s) 



(s''' + 4 /i-y — i^jB s — 3 fi-^C i 



(25) 



Fig. [4]^a) shows the spacing distribution functions 
p^^\s) calculated from the simulation of the point- 
island model and the results of the statistical solution of 
Eq. (19) with 7 = 3,4 and 5. In all cases the agreement 



is good for small values of s and still good for interme- 
diate values of s with 7 = 3 and 4. However, none of 
them describe correctly p^^^ (s) for large values of s (see 
Fig. [4jb)). To improve the model we choose 7 = 4 for 
s < 1.7 and 7 = 3 for s > 1.7, as prompted by Fig. |3] 
This two-regime model, hereafter termed TRM, gives an 
excellent fit with the point-island model even for large 
values of s (see Fig. [4| . 

Fig. reveals the important result that p^^\s) 

does not decay as a Gaussian but instead more like 
exp(— Cs^). This form differs from the ones proposed 
previously in Refs. [UIlT] but is consistent with Eq. (|24|. 

In Fig. [5] we see that Eq. ([24| gives excellent results for 



p'^^^s) even for large and smah values of s. In Fig. ia) 



we show the results of a least-squares fit to calculate the 
extra parameter of Eq. (24). In Fig.jsjb) we use Eqs. ( [25| ) 
to calculate p'^{s). The agreement with the numerical 
data is excellent. 



V. HIGHER SPACING DISTRIBUTION 
FUNCTIONS 

Now we turn to higher-spacing distribution functions. 
Equation (|4| and the weak correlation between gaps sizes 



yields the standard approximate equation for P{s) 



Pis) 



dxp^°\x)p^°\2s-x), 



(26) 



This equation was used satisfactorily in Refs. [TJ [17] to 
calculate P{s) from p'^^^s). Numerical evidence there 
supported the validity of this approximation, which is 
usually called the independent interval approximation 
(IIA). The solutions given by Eqs. a nd (|25]) wih be 
called the analytical model. In Ref. |25j there is addi- 
tional numerical evidence confirming the validity of the 
IIA in this kind of model. In IIA the expression analo- 
gous to Eq. (18) for p'^^^s) is 



«^^+3pW(.) = -<?!(.) + 4,„^(.), 



where we use 



and 



ql{s)= r dxp''^\s-x)pl{s) 
Jo 

q^{s)= r dxp^^\s-x)p^{s). 
Jo 



(27) 



(28) 



(29) 



Equation ( |28[ ) represents the probability density that a 
new nucleation event occurs inside a gap (such as the one 
shown in Fie:. [1]), while Eq. ^ is the analog of p^ (s) 
for this kind of gap. In the IIA q^is) and {s) are 
written as convolution products. The p^^\s) shown in 
Fig. [5] was calculated by using Eq. ([24| and the IIA. We 
find excellent agreement even for small and large values 
of s. There are some small discrepancies (less than 6%) 
near the maximum of p^^\s). 

It can be shown that in IIA, the spacing distribution 
for small values of s can be written as p'^^^s) (x s^^ with 
Qf/c = 3n + 2. In particular, this means that P{s) oc . 
This result agrees with BM's results but differs from the 
exponent reported in Ref. [17^ . 



VI. PAIR CORRELATION FUNCTIONS 

The island-monomer pair correlation function g{r) is 
given by 



iVi(r) 
iVi 



(30) 



where A^i (r) describes the monomer density at a distance 
r from an specified island. From this definition of g{r) it 
is clear that g{r) (x r for small values of r and g{r) 1 
for r ^ 00. In Refs. [ij |26], Ni{r) was calculated by 
using a diffusion equation corresponding to modeling the 
sea of islands as an uniform sink. In this way, they found 



g{r) = 1 - e 



(31) 
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FIG. 6. (Color online) Pair correlation function G{r). 



with ^ the average distance a monomer travels be- 
fore being captured by an island or another monomer. 
The island-island pair correlation function is defined as 
G(r) = N{r) /N, with N{r) the density of islands at a dis- 
tance r from the center of a specific island. In this case 
we have G{r) cx for small values of r; G{r) naturally 
satisfies the same condition as g{r) for large values of r. 
It is well known [27l|28] that, at least in two-dimensional 
systems, G{r) and g{r) satisfy the approximate relation 



G(r) 



(r). 



(32) 



If we extrapolate this observation to our ID case, we have 
G(r) ~ g'^{r). By using Eqs. (31) and (32) we can obtain 



an expression for G(r). These equations predict that g{r) 
and G(r) grow monotonically. However, our numerical 
results, displayed in Fig. |6) show a weak oscillation near 
r = 1. This oscillation has been observed experimen- 
tally in two-dimensional systems [28 - 30 . Both models — 
the one given by the statistical solution of Eq. 
our TRM and the one given by Eq. ([24|) plus 1 



19) with 



]he IIA— 

reproduce correctly this oscillation. The behavior of G(r) 
is shown in Fig.|6] In the analytical model we use Eqs. ([3| 
and (24). The solution for G(r) obtained from the sta- 
tistical solution of Eq. ([l9| is plotted as well. The agree- 
ment between these calculations and the numerical sim- 
ulation is evident. The calculations even reproduce the 
weak oscillations. 



In a test of the validity of Eq. (32), measuring g{r) 



directly from the simulation, we did not find good agree- 
ment with the numerical results. Since G(r) ~ 
arises from the assumption that dN{r)/dt cx A^i(r)*+^, 
apparently the law of mass action is not reliable for the 
whole ring, even though it is valid for the gaps. 



VII. VIABILITY OF GENERALIZED WIGNER 
SURMISE (GWS) FOR EXPERIMENTAL DATA 

A. Comments from preceding analysis 

While the main goal of this paper has been to fully 
characterize the spacing distribution of point islands in 
ID, particularly their tails, an underlying motivation of 
much of our work has been to glean physical information 
from experimental data. For that, the number of real- 
izations are generally less than 1000, often significantly 
less. Hence, the noise in the tails is too large to assess 
the effects discussed above. Fits to data are largely de- 
termined by the range 0.5 < s < 2.5. In that range, the 
agreement between Eq. ([2|) and our numerical results is 
excellent, as seen in Fig. [7| For the ID problem with 
i = 1 treated here, the fitted values of /3 are 1.5 and 4 for 



(0) 



{s) and P{s)^ respectively (for details see Ref. pT]). 
As is especially clear from the inset, the fit in the tails is, 
unsurprisingly, not very good. According to the mean- 
field-like argument in Ref. [17 , the capture zone (CZ) 
distribution can be described by a Fokker-Planck equa- 
tion which has with (3 = 2(i + 1) in ID, as its 
stationary solution. In the associated Langevin equation 



ds 
dt 



K 



(33) 



there is a repulsive force K j3 / s and an attractive force 
proportional to s. Here, is a kinetic coefficient and rj 
arises from the random component of the external pres- 
sure. The latter is responsible for the Gaussian tail of 
P{s) while the repulsive force dominates the behavior of 
P{s) for small values of s. The assumption of Gaussian 
decay underestimates the attractive force which leads to 
a consequent underestimation of the repulsive force. For 
example, in the case oip^^\s) this leads to /3 = 1.5 rather 
than value (3 = 2 given exactly in Eq. ([9| and within error 
in Fig.^b). 

While the GWS does not describe well the tails of 
the observed distributions, it does a considerably bet- 
ter in the more central part of the distribution (typically 
0.5 < s < 2); that part is more germane to analysis of 
experimental data since the data in the tails is often so 
sparse (often below the proper 5 hits per histogram bin) 
that noise makes its use suspect for quantitative analysis. 

We address two questions about the applicability of the 
GWS. First, is there a significant difference in the values 
of P calculated from the whole range of data and from the 
reduced interval 0.5 < 5 < 2? We find /3(an) = 4.229 for 
the entire set of data and /3(sig) = 4.222 for the reduced 
range. Since this difference is insignificant, P can be cal- 
culated from just the data in the central region. Second, 
is the GWS is a reasonable approximation of P{s) in the 
range 0.5 < s < 2, especially in comparison with other 
analytical models? While there is no unique measure of 
"goodness of fit" , we adopt as our quantitative diagnostic 
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FIG. 7. (Color online) The behavior of P(s) and p^^\s) are shown in (a) and (b), respectively. It is clear that the GWS is an 
excellent approximation to these distributions for intermediate values of s. 



the standard Pearson's reduced chi-square statistic [3T] : 



N 



n - 1 ^ 



Pexp(^^.) 



(34) 



where P°^®(s) is the probability of the capture zones with 
size s observed from the numerical data, and P®^p(s) is 
that expected from the analytical model; in the prefactor, 
distinguishing this reduced from x^, A/ is the number 
of bins (typically 239) and n is the number of fitting 
parameters. Another possible metric, e.g., would be the 
area between each model curve and the numerical data, 
i.e., the sum over |P°^'(5^) - P^''p(s^)|. 

In Table I we use to compare fits obtained from four 
different analytical models for the capture zone distribu- 
tion in ID. We consider the GWS with /3 = 4 suggested in 
Ref. [19 (GWS-0), the best fit using the single-parameter 
GWS (GWS-1), the MEM+IIA result given by Eq. ([26]), 
the single-parameter gamma distribution (FD) 



the generalized gamma distribution (GFD) 



A(3S^e- 



(35) 



(36) 



and finally our TRM. Note that with = 2 Eq. ([36) 
reduces to the GWS while for = 1 it becomes TD. 
Clearly, the reduction factor in transforming to in 



Eq. ( 34 ) is here a trivial multiplicative factor since N ^ 
n. However, in addition to computing x^ using the whole 
range of the data, we also find x^ using just the data in 
the interval 0.5 < 5 < 2. In this process, the reduction 
factor is significant in making subsequent comparisons. 

Based on Eq. ([34| the GWS is better than the TD 
and the MEM+IIA but worse than the GTD and the 



TABLE I. Values for the x^(all) and x^(sig) for four different 
analytical models in ID. In the case of the FD we have used 
/3 = a - 1; thus, a ^2/3 + 1 of GWS-1, as noted in Ref. [32]. 





TD 


GWS-0 


GWS-1 


GTD 


MEM+IIA 


TRM 


10^ (all) 


5.155 


0.429 


0.145 


0.038 


1.691 


0.161 


10' x'(sig) 


2.483 


0.353 


0.218 


0.0285 


2.076 


0.162 




1 


2 


2 


1.815 


NA 


NA 




8.612 


4 


4.229 


4.650 


NA 


NA 



TRM. Additionally, the GWS gives a value of (3 similar 
to that of GTD. We recall that in the MEM+IIA model, 
P{s) ^ s^, which is closer to the GTD result than the 
/3 of GWS. With the extra adjustable parameter, the 
GTD gives the smallest x^, but at the price of non-integer 
values for P and u. Not only is the TRM excellent overall, 
it has nearly identical x^ for the full and central ranges 
of 5. It is curious that for GWS-1, 
central range than for the full one. 



X^ is larger for the 



It is noteworthy that the mean-field reasoning of 
Ref. [17 succeeds in ID. This mean- field argument is 
based on knowledge of n(r), which in ID is n{x^y) = 
{2R)-^x{y-x). Similarly, Shi et al.[8 found that in ID 
P4{s) only modestly underestimated the asymptotic peak 
height of the distribution for their two versions of the 
point island model with irreversible growth. The peak 
of P{s) obtained from the TRM is approximately 1.29, 
which is close to the value 1.31 reported in Ref. [8 for 
the limit R ^ oo. 
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B. Extension to 2D 

The maximum entropy method can also be apphed in 
the case of deposition in 2D. We have two normahza- 
tion conditions on the capture zone distribution P{s): 
dsP{s) = 1 and ds s P{s) = 1. We can find an 
additional condition: It is long known that in the ag- 
gregation regime, the densities of monomers and islands 
evolve as A^i ^ and N ^ 0^ ^ respectively [33]. On 
the other hand, the density of monomers, n(r), inside of 
a circular capture zone with radius Vc is given approxi- 
mately by [16] 



n(r) = 0.25 R'^ {rli - r^) + 0.5 R'^ In 



risi 



(37) 



where Visi is the radius of the island and r is the distance 
from the cell center. As before, R = D /F. Then, neglect- 
ing logarithmic corrections, the number of monomers in- 
side a capture zone with radius Vc is proportional to r^, 
i.e., the square of the area of the capture zone. Conse- 
quently A^i - {S'^)N. The analog of Eq. ^ for this 
case is 



(38) 



where jl2 is the second moment of P{s)^ which is constant 
in the aggregation regime. Finally, since P(0) = 0, it is 
reasonable to write 



^MEM 



As^e 



-Cs 



(39) 



where A, B and C are constants. Fig. [Sj shows a com- 
parison between Eq. (39) and numerical simulations of 



island nucleation in 2D. According to Ref. [19 , as well as 
Ref. [18 , Psis), the GWS with /3 = 3 = i + 2 (rather than 
the original mean field result /3 = 2 = i + lin Ref. [17J) 
is a good approximation for P(s); hence, we include it in 
Fig.js] In 2D the GWS evidently approximates P{s) well 
for intermediate values of s but differs substantially from 
the numerical simulations for small and large values of s. 
However, Eq. (39) gives an even better approximation for 
P{s) over the whole range of s (see Fig.|8|. Based on the 
numerical results of Ref. |T8L we use /3 = 4 to calculate 

the Pu^ui^) shown in Fig.js] 

For the GWS, we compare again the values of /3 cal- 
culated from the fits using the whole range of data and 
the restricted range 0.5 < 5 < 2: /3(all) ^ 3.07 and 
/3(sig) ^ 3.04. Again, we evidently can calculate f3 just 
as well using only the data from the restricted range. In 
Table II we quantify with the reduced Pearson's the 
goodness of the fits obtained from some analytical mod- 
els for the capture zone distribution in 2D. In particular, 
we consider the GWS with /3 = 3 suggested in Ref. [19], 
Eq. (39) given by the MEM, and the generalized gamma 



distribution (GFD). For the GFD, we consider not only 
the best fit (yielding f3 ^ 3.86 and u ^ 1.59) but also 
the particular parametrization /3 = 4 and u = 1.5 sug- 
gested by the Evans group ^18| and hence denoted GFE. 



3- 




FIG. 8. (Color online) Capture zone distribution in 2D with 
1 = 1. The GWS describes correctly the behavior of P{s) 
for intermediate values of s. The maximum entropy method 
gives an excellent approximation for P{s) even for large and 
small values of s. 



The one-parameter gamma distribution (FD) is also in- 
cluded. Again we compute x^(all) using the whole range 
of the data and x^(sig) using just the data in the interval 
0.5 < s < 2. Our results are summarized in Table II. 



TABLE IL Values for the x^(all) and x^(sig) for four different 
analytical models in 2D. 





FD 


GWS-0 


GWS-1 


GFE 


GFD 


MEM 


10^ (all) 


3.010 


1.660 


1.726 


0.402 


0.334 


0.518 


10^ x^sig) 


1.722 


0.826 


0.873 


0.381 


0.287 


0.294 




1 


2 


2 


1.5 


1.585 


NA 




6.277 


3 


3.065 


4 


3.860 


4 



Clearly the GWS is inadequate outside the range 0.5 < 
5 < 2, the x^isll) given by the GWS is around five and 
four times as large as the ones given by the GFD and 
GTE{u = 1.5), respectively. In the 2D case, the values 
of x^ of the GWS are closer to the ones given by the 
best fit possible of GFD than it does in ID. However, the 
discrepancies between the GWS and the numerical data 
yield an underestimation of the key parameter p. In fact, 
the difference between the accepted values of f3 and the 
ones given by the GWS in 2D is 25%. The x^(sig) for 
the GWS is around 3 and 2 times as large as the ones 
given by the GFD and GFE, respectively. This suggests 
that the GWS can be used as a first approximation for 
the CZ distribution in the range 0.5 < s < 2. One of the 
main goals of the experiments on epitaxial growth is to 
determine the "critical nucleus" size i from the experi- 
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mental data. Then, in order to estimate i, it is necessary 
to assume an approximate functional form for P{s) and 
then make a fit to find an approximate value for i. Un- 
fortunately, the functional form of P{s) depends on z in a 
non trivial way. The GWS can be used as a first approx- 
imation to determine i from experimental data because 
most of these data are inside the range where the GWS 
is a reasonable approximation. 

We emphasize that the different analytic fitting func- 
tions must be expected to yield different values of /3, 
i.e. the exponent associated with the power-law rise for 
5 < 1. In Table [n] the values of (3 for the two GWS fits 
is considerably smaller (by about 1) than the ^^s from 
Gr and MEM functions, even though all provide decent 
accounting of the central region. The functional forms 
incorporate different underlying assumptions, and only 
the GWS has been related, if approximately, to the crit- 
ical nucleus i. Hence, there is considerable subtlety to 
extracting i from the values of P obtained by fits to ex- 
pressions other than GWS. 

In 2D x^(sig) is invariably smaller than x^(all). How- 
ever, for GFD and GFE, they are relatively close, indi- 
cating the most consistent fit over the whole range, while 
the GWS is significantly better in the central region. The 
GWS is notably better than the FD, often used to ana- 
lyze froths [34 and quantum dot distributions [35] . 

In contrast to ID, the mean-field reasoning of Ref. [17], 
with nucleation probability ex n^(r) (and constant n(r)), 
is not viable in 2D. In 2D the better estimate of n(r) 
given in Eq. (37) comes from solving the diffusion equa- 
tion; however the radial extent of a CZ fiuctuates sig- 
nificantly (as a function of polar angle) around its mean 
value Tc = (S'/tt)^/^. Furthermore, the island is not usu- 
ally at the geometric center of the CZ. Hence, Eq. (37) 



generally does not adequately approximate the local den- 
sity of monomers inside a CZ. In ID newly nucleated CZs 
obtain their "area" from just two existing CZs, while in 
2D several existing CZs often contribute. Most of these 
nucleations occur near the CZ boundaries, which gener- 
ally are not well described as circles. 



VIII. CONCLUSIONS 

BM [1] proposed that p^{y) / FM{y) oc with 7 = 5. 
However, we find that 7 = 3 and 4 give much better 
agreement with our numerical simulation. For example, 
if we choose 7 = 4 from Eqs. ( [12] ) and ( 24 ) it is possible to 
write an analytical expression for the probability density 



Pn^{x,y) in terms of p^^\s). In Fig. [9|we show the re- 
sults of calculating p^^(x, y) numerically compared with 
our analytical model. The contour plots are quite simi- 
lar. However, there is a small difference in the location 
(x, y) of the maximum value of {x, y): In the analyt- 
ical model it is located at approximately (0.7,1.5), while 
in the numerical simulation it is at (0.7,1.3). 

However, neither 7 = 3 nor 7 = 4 describe the statisti- 
cal behavior of the system for large values of s. Further- 



more, the TRM gives an excellent quantitative descrip- 
tion of the statistical behavior of the point-island model 
even for large values of 5, as shown in Fig.[4j This means 
that a complete description of the system requires to take 
into account that 7 is a function of 5, i.e., the probabil- 
ity to choose a gap p^ {y) is different for large and small 
gaps, see Fig. [3) 

We find that the spacing distribution p^^\s) decays 
like ex.p{—Bs^) instead of exp(— ^5^) as in BM or 
exp(— ^5^) as in Ref. [17 . It seems that this is con- 
sequence of the confinement of the diffusive monomers 
within the gaps. An analysis of Eq. ( [T9| ) shows that the 
tail depends exclusively on the probability to choose the 
gap Pn{y)' This suggests that the statistical behavior 
of the system for large values of s is dominated by the 
breakup of the biggest gaps where 7 ~ 3. Meanwhile, 
the behavior of p'^^^ {s) for small values of s depends on 
the probability to choose the position of the nucleation 
inside the gap p^(A). 

Fig. [5] shows that the maximum entropy principle can 
be used to find an excellent analytical approximation for 



p^^\s). Equation ( |24[ ) gives an excellent fit of the numer- 
ical data, even for large and small values of s, in fact bet- 
ter than other approximate expressions [TJ [17] . We find 
that by using the HA and Eq. ( 24 ) , we can easily handle 
expressions for p^^\s) with k > 0. Those expressions ac- 
curately describe the numerical data, even the weak oscil- 
lation in G{r). Apparently, for small values of s the spac- 
ing distribution function satisfies p'^^^s) oc s' 



3fc+2 



We 

show that Eq. ( |18[ ) and ( |24[ ) provide an excellent approx- 
imation for p^ (s)^ as seen in Fig. [s] While ( A) oc (A) 
leads to G(r) oc g'^{r) for small values of r, this relation 
is not satisfied for values of r near the maximum of G(r). 

From our previous results, it is possible to calculate the 
capture number, (jg, of an island with size s in the ag- 
gregation regime. As usual mil2], cr^ can be determined 
from 



1 f dni{x,y) 



dx 



dni{x,y) 



=0+ 



dx 



x=0- 



(40) 



where A^i is the density of monomers and 
{dni{x,y)/dx)\^^Q+[-] represents the derivative of 
ni{x^y) at the right [left] boundary of the island. By 
using the approximation given previously for ni(x,y), 
we found 



y 



CZ 



(41) 



where y^^ is the length of the capture zone. In the point- 
island model the capture numbers naturally do not de- 
pend on the size of the island. This result also implies 
that the distribution of ds has the same form of p^^^(s) 
and not that of p^^\s) as was suggest in Ref. [11 . There 
the authors supposed that the size of the gap at the left 
and the right of an island are equal. A better first ap- 
proximation for the distribution of as can be obtained if 
we calculate p^^^ (5) from p^^\s) through a convolution 
product, as in Fig. p[a). 
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A major result of our analysis, consistent with 
Refs. [lZl[T8i, is that obtaining an appropriate descrip- 
tion of the nucleation mechanism is the crucial ingredient 
to arriving at an excellent approximation for the spacing 
distribution functions of the point-island model. Finally, 
we emphasize that in spite of its mathematical simplic- 
ity, the GWS with the suitable selection of /3 is a good 
approximation for P{s) in ID and 2D. Because of this, 
it is reasonable to use the GWS to analyze experimental 
results for epitaxial growth. 



Advanced Materials (CNAM). 



Appendix A: Behavior of p^*^^ (s) for small and large 
values of s 

An expansion of Eq. ([l4| around s = shows that 
Pn i^) ^ (3Q///^)g^ for small values of s. By using this 



result in Eq. (19), we found p^^\s) « (15//i^)s2. This 



30 as 



result is a consequence of the fact that p^{X) 
A ^ 0. 

From our numerical results, it is reasonable to propose 
p^^\s) ^ As^ e~^^^ for large values of s, where A, 5, a 
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and p are constants. By using this ansatz in Eq. (14) it 
is clear that for large values of 5, ^ p^^\s). Thus, 

Eq. ( 19 ) takes the form 
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(2 + a) 



0, 



^7 



(Al) 



which implies a = —2, B = {p and p = 7. The 
behavior of p^^^ {s) for large values of s is fully determined 
by the probability to choose the gap to fragment. 
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